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Abstract. We show that it is possible to perturb arbitrary vacuum asymptotically flat 
spacetimes to new ones having exactly the same energy and linear momentum, but with 
center of mass and angular momentum equal to any preassigned values measured with 
respect to a fixed affine frame at infinity. This is in contrast to the axisymmetric situation 
where a bound on the angular momentum by the mass has been shown to hold for black 
hole solutions. Our construction involves changing the solution at the linear level in a shell 
near infinity, and perturbing to impose the vacuum constraint equations. The procedure 
involves the perturbation correction of an approximate solution which is given explicitly. 



1. Introduction 

For asymptotically flat spacetimes with appropriate asymptotics there are several con- 
served quantities which can be measured at spatial infinity. These include the total energy 
and linear momentum, as well as the angular momentum and center of mass. When we 
fix an affine frame at infinity the linear and angular momentum as well as the center of 
mass become three vectors. It is natural to ask whether there are any constraints on these 
quantities imposed by the Einstein equations. The positive mass theorem provides one such 
constraint, namely that the energy-momentum vector is a forward pointing timelike vector. 
In particular this says that the magnitude of the linear momentum vector is bounded above 
by the energy. For the Kerr solutions which describe rotating stationary axisymmetric vac- 
uum black holes, it is true that the angular momentum must satisfy such a bound. It has 
been shown over the past several years by S. Dain [11] and P. T. Chrusciel et al. IH O |8] 
that such an inequality is also satisfied by general axisymmetric black hole solutions of the 
Einstein equations. The paper by X. Zhang [17] proves such an inequality under an energy 
condition involving his definition of angular momentum density, but it appears that general 
vacuum data sets do not satisfy this energy condition. 

The main results of this paper show that there are no constraints on the angular momentum 
and center of mass in terms of the energy-momentum vector for general vacuum solutions of 
the Einstein equations. Precisely we fix an affine frame at infinity and we give an effective 
procedure for adding a specified amount of angular momentum to a solution of the vacuum 
Einstein equations, producing a new solution with specified angular momentum but with 
only slightly perturbed energy-momentum vector. We obtain a similar result for the center 
of mass. Then, by considering a family of initial data near the given one, and by doing 
the construction continuously, we obtain a perturbation with arbitrarily specified angular 
momentum and center of mass, while leaving the energy-momentum vector unchanged. One 
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may think of these results as the pure gravity analogue of the addition to a Newtonian 
system of a very small symmetrically placed mass far from an axis whose rotation imposes 
a fixed amount of angular momentum. Similarly one can think of adding a small mass to a 
Newtonian system which when translated far from the center of mass of the original system 
produces a fixed change in the center of mass of the new system. From the point of view of 
the dynamics of the vacuum Einstein equations we expect the angular momentum that we 
add near spatial infinity to be radiated away and to have little effect on the final stationary 
state of the system. We emphasize that our solutions with arbitrarily specified angular 
momentum and center mass are complete manifolds, and we can arbitrarily specify the 
angular momentum of the perturbed data while keeping the energy-momentum and center 
of mass fixed. Without the completeness condition, there are exterior vacuum solutions with 
arbitrary prescribed energy-momentum, angular momentum, and center of mass, such as 
the example of a boosted slice in an exterior Kerr solution computed by Chrusciel-Delay [6] . 
Also, certain A^-body solutions constructed by Chrusciel-Corvino-Isenberg [5] should have 
large angular momentum, but this comes from orbital angular momentum, i.e. c x p for large 
c. In particular, the center of mass is not fixed in their case. 

From a technical point of view the reason it is possible to make these constructions is 
that the angular momentum and center of mass are determined by terms in the expansion 
of the solution which are of lower order than those which determine the energy and linear 
momentum. The idea then is to make perturbations near infinity which affect only the lower 
order terms in the expansion. We do this by explicitly constructing linear perturbations 
supported in a shell near infinity which impose the required change in angular momentum 
(or center of mass), and then by finding a solution of the vacuum constraint equations 
which is sufficiently close to the perturbed system so that the change in angular momentum 
(or center of mass) persists. This can be done in such a way that the energy and linear 
momentum are changed by an arbitrarily small amount. 

Let {M,g,7i) be asymptotically flat in the sense that, outside a compact set, there exists 
an asymptotically fiat coordinate system {x*} so that 

gij{x) = Sij + 0(|x|"^) Tiijix) = 0{\x\~^), 

d''gij{x) = 0{\x\~^-'') for k = 1,2 dn.jix) = 0(1x1"^). 

In addition, we assume that (M, g, vr) satisfies the Regge-Teitelboim condition 

gij{x) - gij{-x) = 0{\x\~^) T^iji^) + 7rij{-x) = 0{\x\~^), 

d\g,,{x) - g,,{-x)) = 0{\x\-^-') for = 1,2 d{7r,,{x) + 7r,,{-x)) = 0{\x\-^). 

The notation / = 0(1x1"°") means that |/| < C|x|~" for a constant C. We remark that our 
construction works for data [g, it) with weaker assumptions on the decay rates. For simplicity 
of notation, we assume the decay rates above and do not consider here the question of optimal 
decay conditions. 
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Let -E, C, P, J denote the energy, center of mass, linear momentum, and angular momen- 
tum of {g, vr). They are defined as limits of integrals over Euclidean spheres 




where dao is the area measure of the Euclidean sphere {\x\ = p} and Yi = -£j x x (cross 
product) for z = 1, 2, 3 are the rotation vector fields. Denote by E, C, P, J the energy, center 
of mass, linear momentum, and angular momentum of (^, vf). 

We now give precise statements of the main theorems, where the definition of the weighted 
Sobolev spaces W'l'g is provided in Section [3l and we assume p > 3 and q G (1/2, 1). In our 
construction, we fix an affine frame near infinity and measure all asymptotic quantities rela- 
tive to this frame; in fact, we may fix an asymptotically flat coordinate system throughout. 

Theorem 1. Let (5',7r) be a nontrivial vacuum asymptotically flat initial data set with g = 
outside a compact set and 7r(a;) + 7t{—x) = 0{\x\~^~'^'^) . Given a G and e > 0, 
there exists a vacuum asymptotically flat initial data set {g,7f) such that {g,n) is within the 
e-neighhorhood of {g, tt) in Wlg x W^f_g and 

\E-E\<e, |C-C|<e, |P-P|<e, 

and 

(1.1) |J-J-a|<e. 

For the center of mass we prove the following. 

Theorem 2. Let {g,iT) be a nontrivial vacuum asymptotically flat initial data set with g = 
outside a compact set and tt{x) + 7r(— x) = 0{\x\~^^'^'^) . Given 7 G and e > 0, 
there exists a vacuum asymptotically flat initial data set i^iTf) such that (^, tt) is within the 
e-neighborhood of {g, vr) in Wl'^ X W^f_g and 

\E-E\<e, |J-J|<e, |P-P|<e 

and 

(1.2) |C-C-7|<e. 

By combining these two results we can change both the center of mass and angular mo- 
mentum so that they are arbitrarily close to specified values while leaving the energy and 
linear momentum essentially unchanged. The condition that g = v^6 can be removed by a 
density theorem. Moreover, given a vacuum initial data set, there is a small perturbation 
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with arbitrary specified angular momentum and center of mass and witli the same mass and 
linear momentum. 

Theorem 3. Let {g, vr) be a nontrivial vacuum initial data set satisfying the Regge-Teitelboim 
condition. Given any constant vectors ao;7o ^ there exists a vacuum initial data set 
{g, vf ) within a small neighborhood of {g, vr) in Wl'!; X W_f^g and 

E = E, P = P, 

and 

J = J + tto, C = C + 7o. 

In Section 2 we give an explicit construction of solutions of the linearized constraint 
equations which satisfy a certain moment condition. Sections 3, 4, and 5 are devoted to the 
proofs of the main theorems. We remark that the constructions of Section 2 are explicit, 
while the method of solving the exact constraint equations from the approximate solution 
involves constructing a small solution of an elliptic system with leading order term the 
diagonal Laplace equation. It should be possible to numerically approximate the resulting 
solutions to a high degree of accuracy. 

2. Compactly supported solutions of the linearized constraints 
Recall that the vacuum constraint equations for initial data {g, vr) may be written 

R{g) + i (Tr.vr)' - |vr|2 = 0, div,(vr) = 0, 

where Hij = Kij — TTg{K)gij is the momentum tensor. In general, we consider the constraint 
map $ defined by 

$(^7, vr) = {R{g) + ^ (Tr.vr)^ - |vr|2, div,(vr)). 
The vacuum constraint equations, linearized at the trivial data (5,0), become 

La := ^{(Jij,ij - aujj) = 0, div(r) = 0, 

for symmetric (0, 2) tensors (a, r). In this section we will construct solutions of the linearized 
constraint equations which are compactly supported in the shell Ai = {x : 1 < |x| < 2} 
contained in and which have certain specified moment conditions with respect to rotation 
vector fields. We use the Einstein summation convention and sum over repeated indices; 
though, sometimes we employ summation symbols for clarity. 

We write the Euclidean metric on in spherical coordinates dr^ + r^gabdx°'dx'' where g is 
the standard round metric on S^. The coordinates are labeled by r = x° and 6,4> = 
The ranges for the indices are i,j,k,l,... = 0,1,2, and a,b,c,d,e, ... = 1,2. If a, (3 are 
one-forms, we define the symmetric product a Q (3 to he the symmetric (0, 2) tensor whose 
components are 

(a f3)ij = ^{aif3j + ajf3i). 
We first impose the following ansatz for our solutions. 
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Lemma 2.1. Suppose that q and Q are functions of r, a is a one-form on S"^, and f is a 
trace-free symmetric (0, 2) tensor on S"^ . Then r = 2qa Qdr -\- Qf is a trace-free symmetric 
(0,2) tensor on Ai. The condition div(r) = becomes 

diva = and (r'^qYa + Q(divf) = 0, 

where div is the divergence operator of 5^ on tensors. 

Suppose that p and P are functions of r, fj is a one-form on S^, and a is a trace- free 
symmetric (0,2) tensor on S"^. Then a = 2pfj Q dr -\- Pa is a trace-free symmetric (0,2) 
tensor on Ai, and La = if 

2r(rp)'divr/ + P(divdiv5") = 0. 

This lemma follows directly from the following two computational lemmas. For a coordi- 
nate system on M^, denote by gijdx^dx^ the Euclidean metric. 

Lemma 2.2. Let h he any symmetric (0, 2) tensor on M? , then 

Proof. Let V be any vector field. We have V\h\)-k = {y''h^)\k — h^V-l- Now 

Therefore V\h^^),k = V'^^kiy^h'^) - h'^TliV^ and then 

Lastly, we plug in the formula for T^j. □ 
Now we apply this to the spherical coordinates 

Qijdx'dx^ = dr'^ + r'^dO'^ + sin^ Odcj)^ = dr'^ + r'^gahdx°'dx^ 
and write h = ho^dr'^ + 2hoadx'^dr + habdx'^dx'^ and a = aodr + aadx°-. 
Lemma 2.3. Let h he a symmetric (0,2) tensor and a a one-form on M? , then 

and 

We shall take the following ansatz on h that Hqq = and g"'''hab = 0. We also assume that 
= p{r)fja and hab = P{r)hab for a one-form fj and a symmetric (0, 2) tensor h on 5*^; that 
is, we have 

h = 2p{r)fi Q dr + P{r)h. 
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We then have 

divh = {r"'^p) {div fi)dr + r'''^{r'^p)'fi + r~^Pdivh 

and ^ ^ ^ 

divdiv/i = 2r^'^(rp)'(divr/) + r^^P (divdiv/i). 

Hence, Lemma [2.11 follows directly from the above two identities. 
In particular, if we consider the tensor r of the following expression 

(2.1) r = -2gdi^(f ) rfr + (r2g)'f . 
Then r satisfies divr = on if 

(2.2) di^di^f = on 
Also, we consider a as follows: 

(2.3) a = — 2pdiv5" Q dr + 2r{rp)'a. 

Then a satisfies La = on Ai for any trace-free symmetric (0, 2) tensor 5" on S*^. 

In the following, we show that there are nontrivial solutions to (12. 2p . We can start from 
any one-form fj on S"^ and construct a trace-free symmetric (0, 2) tensor 5*7] on S*^ by defining 

{5*Vl)ab = ^{Va;b + Vb-a " g''%;d9ab)- 

We need the following computational result. 
Lemma 2.4. 

div5*f/ = --{dd* + d*d)fi + f]. 
In particular, if u is a function on S'^ , then 

diyS*du = l:d{Au + 2m), 

and 

disr5*{*du) = 2 * d{^'u + 2m). 

Proof. We recall that for one-forms on S*^, we have d* = — * d*, and for a function u, 
d*du = —Am. We assume fj = du and compute in an orthonormal frame: 

(div(5*dM))i = ^(M1;11 - M2;2l) + Ml;22- 

On the other hand, {^{dd* + d*d)du)i = — |(An)i = — |(ni;ii + U2-2i)- Therefore, 

{diyS*du)i + (-(c?(i* + d*d)du)i = ui-22 — M2;2i- 

At last, we use the commutation formula Ub-ab = Ub-ba + Ua ioi a ^ b on S^. Other formulae 
can be checked similarly. □ 

The next lemma shows how solutions of Equation (12. 2p can be constructed: 

Lemma 2.5. Suppose that a is any one-form on S^. Then f = 6*a is a trace-free symmetric 
(0, 2) tensor on S'^. Moreover, if a = *du for a function u on S"^, then f satisfies divdivf = 0. 
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Proof. This follows from the formula of Lemma 12.41 

divf = * d{Au + 2m), 

together with the fact that (d*)^ = 0. □ 

We need the following computational result. 

Lemma 2.6. Suppose that t = 2qa Q dr + Qf, a = 2pfj Q dr + Pa, and Y = Y"-^ is 
tangent to S"^. Then 

lr,j,iY^ + tuY;^ 

= pqr-\a,,aY'' + aaK^^c^'' + \PQr-\hc;aY'' + r^feX^ + facY-^Yg'^'a^.-g'^ 

where Ob-a, Tbc-a, o,nd Y.^ denote covariant derivatives of a, f, and Y with respect to the 
standard metric gat on S"^. 

Proof. Direct computation. □ 

We are finally in a position to prove the main results of this section. 

Theorem 2.1. Given any A = (Ai,A2,A3), there exist symmetric (0,2) tensors a,T E 
C^{Ai) with (j{x) = a{—x) and r(x) = t(— x) satisfying 

La = 

(2.4) J]r,,„ = 0, /orj = 1,2,3, 



so that 



(2.5) 



o"*-' dx = Afc 



for k = 1,2,3 where Yk = x x (cross product). 

Remark. By a direct computation, the integrand in (12. 5p equals |(iiy^r)jj(T*-', where -Cy^r 
is the Lie derivative of t along Y^ on . Thus, if t is axisymmetric with respect to Y^, i.e. 
S^Yk^ = 0, then (12. 5p is always zero. 

Proof. We first show how to make the integral on the left of (12.51) for k = 1 nonzero. To 
simplify notation for this purpose we denote Yi hj Y. We choose r of the form (12.11) for 
some f = S* * du, where u is an even function on S"^. By Lemma 12.11 and Lemma 12. 5[ r 
satisfies Equation (12. 4 p and has the desired symmetry r(x) = t{—x). We take a of the form 
(12. 3 p for a symmetric (0,2) tensor a to be determined later. By Lemma \2.6\ the integral in 
question can be written as 

" pqr-^dr) I (a^^^F" + aaY.l)f,c9''dS^ 
+ ( / IPQr-'dr) [ {nc;aY^ + f,,K^ + facY.tyg^'ader'dS^ 
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for a = — divr, t] = —diva, P = 2r{rp)', and Q = (r'^q)'. We can choose p and q to be any 
compactly supported functions on the interval (1, 2) so that (J^ pqr~'^dr) and {J^ ^PQr~^dr) 
are arbitrary. It suffices to choose a to make the following integral nonzero 

To achieve this, we take 

and define a by Equation 02.31) . Since Y is invariant under x H- — x, a defined in this way 
has the desired symmetry and satisfies La = 0. It is not hard to check that a = -Cyf, the 
Lie derivative of f with respect to Y on S"^. Thus, we can take an even function u (e. g. the 
restriction of any homogeneous polynomial of even degree) so that fbc-^oY'^ + fabY.'^ + facY.l 
is nonzero. 

To achieve the desired conclusion, we consider the linear functional T(o-,t-) (^^0 given by 
the left-hand side of (12. 5p with vector field Y = v ^ x. Since this is a nonzero linear 
functional we may choose a positively oriented orthonormal basis {61,62,63} so that the 
vector (T(o-T-)(6i), T(o-T-)(62), T(o-T-)(63)) is proportional to A, and after multiplication of r by 
a constant we may assume the vector is equal to A. It follows that there is a rotation R 
of so that T(^,^)(i?(^)) = A^ for k = 1,2,3. It follows that (E^D holds for the pair 
{{R~^)* (a) , {R~^)*{t)) since we clearly have T(^s''a,s*T) = ^(o-,t) ° for any rotation S. 

□ 

We will need a corresponding result which will be used to specify the center of mass. This 
involves the construction of solutions of La = satisfying a moment condition. 

Theorem 2.2. Given any /3 = /32! /^s) G M.^, there exist a trace-free and divergence-free 
symmetric (0,2) tensor a G C^{Ai) satisfying La = so that 




forp= 1,2,3. 

Proof. We first show how to make the integral on the left of (12.61) nonzero for p = 1. By 
starting with a nonzero function u supported in the first octant, we find from Lemma [2.51 
and (12.11) a nonzero trace-free symmetric (0,2) tensor a in C^{Ai) satisfying div(o") = 0. 
Then, in particular La = 0. Because a is supported in the first octant, this implies 

/ x^y^iaij^kf dx > 0. 
We then let T^- be the linear functional on M"^ given by 




Since is nonzero, there is a positively oriented orthonormal basis {61,62,63} for which 
the vector defined by (To-(6i), To-(62), To-(63)) is proportional to /3. Replacing a with a scalar 
multiple we may assume that To-(6p) = /3p for p = 1,2,3. Thus there is a rotation R 
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with -R(gf^) = Cp SO that T„{R{-^)) = /3p for p = 1,2,3. Since for a rotation S we have 
Ts*a = o it follows that {R~^)*a satisfies the required condition (12.61) . □ 

3. Specifying the angular momentum 

We first state a general analytic result which constructs new initial data sets from given 
ones. Let La := Ylii ji^ij,i3~^ii,j3) be the linearized (at the Euclidean metric) scalar curvature 
map which goes from smooth symmetric (0, 2) tensors to smooth functions. Denote by W^'^ 
the weighted Sobolev spaces defined as follows. We say / G Wtq if 



where a is a multi- index and p is a continuous function with p = \x\ on the region where the 
asymptotically fiat coordinate system {x*} is defined. When p = oo, 

||/||^.,o.= ^e..sup|D"/|pH+'^. 

We assume /c = 1 or 2, g G (1/2, 1), and p > 3. By our assumption [g, tc) G W'^^ x W]_f_q. 

In the following, the notation / = 0(r~") means that |/| < Cr~°-, \df \ < Cr~°'~^ for 
some constant C independent of k and the analogous conditions on successive derivatives 
as needed. We remark that the following results would hold similarly for any p > 3/2, but 
the "O" -notation below would mean the decay in weighted Sobolev norms. Here we assume 
p > 3 because the weighted Sobolev norms can be replaced by pointwise estimates using the 
Sobolev imbedding theorem. 

Let r := Denote the shell hj Ak = {k < r < 2k} and denote by C^{Ak) the set of 
functions or tensors which are compactly supported in Ak- 

Proposition 3.1. Let the symmetric (0,2) tensors a,T & C^{Ai) satisfy the linearized 
constraint equations 

La = 
YTij,i = for j = 1,2,3. 

i 

Define a^,T^ G C^{Ak) by 



a 



k ^a{x/k), and = k "^r^x/k). 



Then given any vacuum initial data set {M,g,Ti) with decay rate g — 6 = 0{r~^), vr = 0(r~^), 
and any fixed q G (1/2, 1), there exists a sequence of vacuum initial data sets {(^*^,vf'^)} so 
that for k large and outside a fixed compact set (independent of k), 



(3.1) g^^=[i + ^]g^^+a^^+Oir 



(3.2) 4 = TT,, + r^. + 1 



-B^Xj - B^Xi + BfxiSij 



0{r 
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where A'' and (-Bf , -B|, -B|) are constants. The initial data sets {(^^',vf^)} are small pertur- 
bations of {g, tt) in the weighted Sobolev spaces; in fact, 

— ^'^||^2,p — )■ 0, ll^r — Tf'^lliA/i.p —^0, as k oo. 

'''^ -q -1-q 

Moreover, 

(3.3) Ek E and P as k ^ oo. 

Proof. Let g^ = g + cr^ and vr^ = tt + r'^. Then {g'^^ii^) satisfies the constraint equations 
(^{g^.Ti^) = (0,0) everywhere in M\Ak and ^{g'',Tf^) = {0{r~^),0{r~'^)) in A^. We denote 



for any vector field X and metric g. By the proof of [TOl Theorem 1], there exist (u^, X^) on 
M and {h^,w^) with compact supports (uniformly in k) such that 



satisfy $(^'^,7f^) = for all k large, and 

(3.4) \\u'' -l\\^2,p ^0, ||X*^||j^2.p ^ 0, ||/i^||^2,p 0, ^ 0, 



-q -q ~q -q 

as A; — 7- oo. The constraint equations imply An'^ = 0(r~^~^^) and A(X^')j = 0(r~^~^'?). It 
follows that 

M*^ = l + 4^ + 0(r-2'?), and CX% = ^ + 0(r-^'i). 

Therefore, (13.11) and (13.21) follow, and the convergence of (13.31) can be derived as in ^U\. □ 

For the rest of the section, we consider the special case of Proposition 13.11 when g = 
outside a compact set so we have = (l + ^) Sij+0{r~'^). We also assume that vr = 0(r^^), 
7r(x) + 7r(— x) = 0{r~^~'^'^). That any initial data can be approximated by such data follows 
from [To]. Note that the asymptotic oddness condition is the Regge-Teitelboim condition 
fl6\ [T] required for the existence of J. Therefore, near infinity, {'g'',w'') satisfy the conditions 

(3.5) 4- = {u'rig., + 4) = (l + 5., + 4 + 0{r~'^') 



(3.6) = TT,, + + ^ 



0(r 



Proof of TheoremUl Let {g, vr) be a vacuum initial data set satisfying the above conditions. 
Choose a and r satisfying the assumptions in Theorem 12.11 with A = —8na. There exist 
vacuum initial data sets {g'^, If^) satisfying (13. 5p and (13. 6 p by Proposition 13. H and \Ek—E\ < e 
and |Pfc — P| < e. It remains to prove the desired properties of the center of mass and angular 
momentum. In the following, we suppress the superscript k of 'g^ and If^ whenever it is clear 
from the context. 
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We assume that a = (ai,a2,Q;3). We denote by Yp the Euchdean Kilhng vector field 
= X X for p = 1, 2, 3. Because of the asymptotics of {g,7r), 

EJp = ^ hm / Y^w.^iYpY^ dao 
= ^ hm / WijiYpYu^ da-g, 



where v and da-g are respectively the outward unit normal vector and the area measure of 
{r = p} with respect to "g. To shorten notation we do the following estimates using Y in 
place of Yp. By the divergence theorem, assuming po <^ k and 2k < p, 

/ 7rijY'uUag= / g^^(WiiY'),jdYolg+ / WijY'uUag. 

J r=p PO<'"<P J i'=po 

Because r'^ vanishes on {r = po}; 

I 7r,jY'u^dag= I Y^iu^fi-K + CgX^),jY'uUag. 

Jr=po Jr=po jj- 

Although the absolute value of the last integrand is 0(po ^) which would not have a finite 
limit, the asymptotic symmetry conditions say that vr is odd and the term CgX.^ is also 
asymptotically odd, and hence the leading order term of the integrand is odd causing the 
limit to be finite. In fact, we have for k large enough, 



Jr=po ij 

= I $^vri,rV^ da^+ I J2{CgX.%jY'iy^ da-g 

Jr=po j j J r=pn „• „■ 



(3.7) 

(3.8) + I Y}^u^^^ - + /:§X^),,rV^- dag 

Jr=po 

= 87rE3p + Oipl'"'). 

Clearly, the first integral in (13. 7p is 8TrE3p + 0{pq^). For the second integral in (13. 7p . we use 
(13. 4p and choose k large so that l^gX'^l is small, say less than pg^. The integral in (13.81) is 
0(pq^^'') by (13. 5p . (13. 6p . and the asymptotic symmetry of tt. 

To estimate the interior integral, by the constraint equation 'g^^Tfu-j = and the condition 
that y is a Euclidean Killing vector field, 

(3.9) 9'^{^iiY% = it' - 5'^>,,yi + g'^7faY%. 

By (13. 5p . (13. 6p . and a^{x) = a^{—x), the integral of the first term on the right-hand side is 

J Pn<r<p JAi. ■ ■ , 



'0 

- [ 5^a,,r^;1^5rfa; + 0(pJ-'^) 
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where we use 7i{x) + -7r(— x) = 0(r ^ ^'') and = a''{—x) to estimate the error terms. 

For example, for some of the error terms, 



2EA^ 

Po<r<p r 



i ^ 1 J pa 



J2 '^ki^'i^J ^^OV = 

PO<r<P ij^i 

To estimate the integral of the second term on the right hand-side of fl3.9p . we again use 
(I3.5p . fl3.6p and asymptotic symmetry to derive the first equality. 

= -\l 44,^' dx + oipi'") = -U J2 dx + o{pi-'% 

where in the second-to-last identity, we integrate by parts and use the fact that and 
vanish on the boundary and that Y is divergence-free. 
Combining the above identities, we derive 



PO<r<p 



aij dx + 0{p, 



J 



Then we choose po large so that the error term is less than e. For k ^ po large enough, we 
prove that {g^^n^) satisfies (11. ip from Theorem 12. II applied with A = — SttcJ. 

We show that the center of mass of (^, vf) remains almost unchanged during the process. 
For p = 1, 2, 3 we have the components of C defined by 



= lim 



r=p 



X' 



1,3 



X 



X^ 



dan. 



By the divergence theorem. 



IGnE = lim 



/ ^^'Y^(9ij,ij 9ii,jj)dx 

J P0<r<p jj 



(3.10) 



Jpo<r<p 



r=po 



Z. 2^ \ 9ip~ 9ii~~ 



daQ. 



Because J2i j{9ij,ij~9ii,jj) is the leading order term of the scalar curvature, it can be replaced 
by the lower order terms such as |7fp and (Dg)^ using the constraint equations. Then by 
(13. 5p . (13. 6p . and the symmetry cr'^(x) = a'^(— x), the interior integral above is 0(po~^''). 
Similarly, we have 



(3.11) IGttECp 



r=po 



^ 2^y9ij,i 9ii,j) ^ 2^ [ 9ip ^ 9ii ^ 



dao + 0{p^'). 
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Moreover, because g = {u^Yg and u^' is close to 1 on {r = po} for k large enough, the 
difference of the boundary integrals on {r = po} is 0{pq^). Therefore, for a fixed po large 
and for k ^ pq large enough, we have 

|C- C| < e. 

□ 



4. Specifying the center of mass 

As in the previous section we assume that g = v'^S outside a compact set and tt = 0(r^^) 
and 7r{x) + 7i{—x) = 0(r~^~^'^). We apply Proposition 13.11 with r = and a chosen by 
Theorem 12.21 to be a solution of La = satisfying the moment condition (12. 6p . 

If. 

Proof of Theorem\^ By Proposition 13. \E — E\ < e for k large. Since r = 0, the proof 
that the angular momentum satisfies |J — J| < e follows as in the previous section. 

To estimate the change in the center of mass we use fl3.10p . (13. lip , and the argument 
following them. For p = 1,2,3 we have 

167r (EC' - EC^) = hm [ x^ - 9u,,) dx + 0{p-,'). 



Because a^{x) ^ o"^'(— x), the interior term above is not of lower order. Let R denote 
the scalar curvature of 'g. Then from the constraint equations, we have R = 0(r~^) and 
R{x)-R{-x) = 0(r^3-2g)^ gy [Y3| Lemma 3.5], 



X' 



po<r<p 



dx 



/ x^Rdx- I xPSgdx + 0{po^ 

J po<T<p Jp()<r<p 



pO<r<p 



where 



^g-~ - ^ii){'^9ij,ij - Qiijj - 9jj,ii) 



>r-r _ _ _ _ 3_ _ 1_ _ 1- - 1 

/ , [ 9ji,j9u,i + 9ji,j9ii,i + 'i9ij,i9ij,i '79jj,i9ii,i 7^9ij^i9ii,j\ ■ 



Recall that a'^ G C^{Ak) is trace-free and divergence-free. By (13. 5p . the above integral is 
equal to the following, up to an error term of order 0(pq~^^), 

'2E + A''^ 



k k 



f2E + 
\ r 



L r'Y. ( 

i,j,i I ^ 



2E + A'' 



dx 



dx 



2E + A^ 



dx 



X' 



),i dx= I xP^j{aij,ifdx 
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where in the last hne we use integration by parts. We may then apply Theorem 12.21 with 
P = 647ri?7 to obtain the required condition (II. 2p on the center of mass. □ 

5. Proof of Theorem [3] 

In Theorem [1] and Theorem [2] we assumed that g = v^6 outside a compact set and vr = 
0(r~^), 7t{x) + 7t{—x) = 0(r~^~^^). Using a density theorem [12], we prove below that the 
condition can be replaced by the weaker Regge-Teitelboim condition. 

In this section, we fix the constant p > 3 and the constant g e (1/2, 1). 

Theorem 5.1. Let {g,TT) be a nontrivial vacuum initial data set satisfying the Regge- 
Teitelboim condition. Given a, 7 G and given e > 0, there exists a vacuum initial data 
set {g, vf) with \\g — 5'|Itia2,p < e, ||7f — ttHtx^i.p < e, so that 

—q * —1—q 

(5.1) \E-E\<e, |P-P|<e, 
and 

(5.2) |J-J-a|<e, |C-C-7|<e. 

Proof. By the density theorem in [12], given the vacuum initial data set {g,7T) satisfying 
the Regge-Teitelboim condition and any e > 0, there exists a vacuum initial data [g, n) 
with g = v'^S outside a compact set and 7f(x) = 0(r~^), 7r(x) + 7r(— x) = 0{r~^~'^'^), so that 
Wq " dWw^-'^ — ^' 11^ ~ ^llw^'f — ^- Moreover, 

\E-E\< -, IP - PI < -, |J - J| < -, |C - CI < -. 

I I - 3' I ' - 3' ' ' - 3' ' ' - 3 

Then we apply the construction in the proof of Theorem[T]to {g, it). Given a = {ai, a2, as) G 
M^, there exist symmetric (0,2) tensors a, r G C^{Ai) satisfying 



(5.3) 



a^^ dx = —Snap, 



where rotation vector fields Yp = ^ x x. Let a''{x) = k ^a{x/k) and t''{x) = k '^a{x/k). 
By Proposition 13. H there exists a large integer k so that [g^.,Ti^): 



satisfies the vacuum constraint equations, where {u'^, X^') and {h'', w^) arise from solving the 
linearized constraint equations. They satisfy the decay condition (13. 4p . and the [h^, w^) have 
compact support. Then, from Theorem [1], we have 

\E-E\<-, |P-P|<-, |C-C|<-, 



and 

e 



3' ' ' - 3' ' ' - 3' 

J- J-a| < 



3 
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\E-E\< -, 


IP-PI < -, 


|J- J| < 


1 1-3, 


1 1 - 3> 





We suppress the superscript k of {g^, tt^) in the following. Clearly (gf, vr) satisfies the condition 
in Theorem [21 namely g = {uv)'^6 outside a compact set and vr = 0(r~^), 7r(x) + 7r(— x) = 
Q^^-i-2qy j^g^ 0" be a symmetric (0,2) tensor satisfying 

(5.4) [ xP^(a,,-fc)2rfx = 647rE7p. 

Let a\x) = l^^a{x/l). By the construction in the proof of Theorem[2l there exists an integer 
I ^ k so that 

~g^ = (uY(g + a') + h^, 

^' = (^f (7r + £(^+^.)X')+u;' 

satisfies the vacuum constraint equations, where ('u',X') and {h\w^) are from solving the 
linearized constraint equations as above. Moreover, 

e 

3' 
and 

Then (15. ip and (15. 2p follow by combining the above inequalities. □ 

We can further perturb (^, vr) so that the energy-momentum vector equals to that of ((?, tt), 
while changing the angular momentun and center of mass by only a small amount. 

Proposition 5.1. Let {g-iH) he a nontrivial vacuum initial data set satisfying the Regge- 
Teitelboim condition. Given a, 7 G M'^ and e > 0, there exists a vacuum initial data {g, vf ) 
satisfying \\g — g\\u^'2,p < e and ||7f — ttHtj^i,? < e so that E = E and P = P and 

I J — J — a| < e, and |C — C — 7I < e. 

Proof. Let {g, vr) be the initial data constructed in Theorem 15. 1[ Let {g, vf ) be the vacuum 
initial data from scaling g = X^g and vf = Avf , where the constant A is a positive constant and 
= (E'^— |Pp)/(£'^ — |Pp). Then, by straightforward computations, we have (£',P, J,C) = 
A(E,P, J,C) and then 

Notice that by (El]), 



E^ - |P|2 = E' 



{E' - |Pn - {E' - |Pn < 2e{E + |P|) + e\ 

Therefore, since -E > |P| by the positive mass theorem, we divide the above inequality by 
E^ - |P|2. Then 

A-2 - 1 < - — — + 



E-\P\ E^ 



12 • 



Because E and P are fixed, we can choose k, I in the proof of Theorem 15.11 large enough so 
that A is close to 1. Therefore, {E, P, J, C) is close to {E, P, J, C) and hence to {E, P, J, C). 
Because E is close to E, we then boost the data {g, vf ) by a small angle so that E = E, 
and then |P| = |P| (the existence of such boosted slice is proven in [5]). By rotating 
the asymptotically flat coordinates, we can make P = P. Also, notice that the angular 
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momentum and center of mass only change a small amount after these transformations. 
(The transformation formulas of these quantities under the Poincare transformations can be 
found in, for example, [6l Appendix E].) □ 

To prove Theorem [3l we need the following degree argument. 

Lemma 5.2. Fix the constant a > 0. Let Baizo) C M" denote the closed ball centered at Zq 
with radius a. Let f : Ba{zQ) — > M" be a continuous map satisfying, for any z G Ba{zQ), 

\f{z)-z\<a. 

Then /^^(^q) is non-empty. More precisely, either f{z) = Zq for some z G dBa{zo) or the 
degree of f at zq is one. 

Proof. By scaling, we only need to prove the case when a = 1. We define the continuous 
homotopy between / and the identity map for < t < 1: 

h{z,t) = {l-t)z + tf{z). 

For a boundary point z G OBi^zq), 

{h{z, t) - Zo) ■ {Z - Zo) = [{Z - Zo) + t{f{z) - Z)] ■ {Z - Zq) 

> \z — ~ 't\f{z) — z\\z — Zo\ > 1 — t. 

Then either f{z) = Zq for some z G dBi{zQ) or h{z,t) ^ Zq for all < t < 1 and for all 
z G dBi{zo). In particular, the latter case implies that zq stays in the range of h{-,t) for all 
t G [0,1]. Therefore, f~^{zo) is non-empty. □ 

Proof of Theorem\^ Denote the given constant vector ((JojTo) by zq . We may without 
loss of generality prove only for the case |ao| ^ and 1 70 1 7^ 0, for if ao (or 70) is the zero 
vector, we apply the theorem twice to a non-zero constant vector v and then to —v. We 
define the map / : 5,(^o) C ^ by 

/(a,7) = (J-J,C-C), 

where J and C are the angular momentum and center of mass of ((?, 7f ) constructed in 
Proposition 15. 1[ By the construction, 

\f{z)-z\<e. 

Once we verify that / is continuous, we apply Lemma [5l2] to obtain f{d,^) = (ao,7o) for 
some (a, 7) and complete the proof of Theorem [31 

Claim. The map f is continuous. 

Proof. Given non-zero vectors (aoiTo); let (cr, r, a) be the symmetric (0,2) tensors and k,l 
be the integers in the proof of Theorem 15. II for {do, 70). (Notice that k, I may be chosen large 
depending only on do,^o, and e.) 

Assume that (a, 7) is another pair of constant vectors. We fix the symmetric (0, 2) tensors 
a, T, a and the annular shells determined by k and /. We choose the symmetric (0, 2) tensors 
for a, 7 from a, r, a and apply the construction in Theorem 15.11 over the annular shells 
determined by k and /. 
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That the choice of the symmetric (0, 2) tensors depends continuously on a. can be seen 
as follows: Let i?i : -> be a rotation so that i?i(ao) is parallel to a. The symmetric 
(0, 2) tensors defined by 




satisfy the corresponding condition fl5.3p for a on the right-hand side. It is easy to check 
that at each point 

' ,^ ^ , / 1 , , 3,^ 

<\oi- aol 7TT^|c^l + ol^^ 
\2|ao| 2 




Similar estimate can be derived for the other tensor r. For 7 and the integral (15.41) . we can 
also choose the tensor corresponding to a in the same fashion. It is straightforward to check 
that the rest of the construction is continuous, and hence / is continuous. 

□ 

Corollary 5.3. Given any constant vector (i?,P, J,C) G M}^ with E > \P\, there exists a 
smooth and complete asymptotically flat vacuum initial data set whose energy, linear mo- 
mentum, angular momentum, and center of mass are the corresponding components of this 
constant vector. 

Proof. By Theorem [31 it suffices to show that there exists a vacuum initial data with the 
specified E and P. By the results of the global existence of the Cauchy problem (see 
[21 El US])) given a strongly asymptotically fiat vacuum initial data set ((?, vr) close to the 
fiat data, there exist future and past complete vacuum developments. In particular, we can 
boost the slice in spacetime and then rotate the coordinates so that the energy-momentum 
vector of ((?, vr) is parallel to the given vector {E,P). Then, by scaling the data, we obtain 
a vacuum initial data with the desired energy-momentum vector [E, P). □ 



References 

[1] Beig, R. and O Murchadha, N. The Poincare group as the symmetry group of canonical general 
relativity, Ann. Physics 174 (1987), no. 2, 463-498. 

[2] Christodoulou, D. and Klainerman, S., The global nonlinear stability of the Minkowski space, Prince- 
ton Mathematical Series, 41. Princeton University Press. 

[3] Christodoulou, D. and O Murchadha, N., The boost problem in general relativity, Comm. Math. 
Phys. 80 (1981), no. 2, 271-300. 

[4] Chruscicl, P. T., Mass and angular-momentum inequalities for axi-symmetric initial data sets. I. 
Positivity of mass, Ann. Physics 323 (2008), no. 10, 2566-2590. 

[5] Chrusciel, P. T., Corvino, J., and Isenberg, J., Construction of N-body initial data sets in general 
relativity, larXiv: 1004.13551 

[6] Chrusciel, P. T. and Delay, E., On mapping properties of the general relativistic constraints operator 
in weighted function spaces, with applications, Mem. Soc. Math. Fr. (N.S.) No. 94 (2003). 

[7] Chrusciel, P. T., Li, Y., and Weinstein, C, Mass and angular-momentum inequalities for axi- 
symmetric initial data sets. II. Angular momentum, Ann. Physics 323 (2008), no. 10, 2591-2613. 

[8] Chrusciel, P. T. and Costa, J., Mass angular-momentum and charge inequalities for axisymmetric 
initial data. Class. Quantum Grav. 26 (2009), no. 23, 235013, 7 pp. 

[9] Corvino, J., Scalar curvature deformation and a gluing construction for the Einstein constraint 
equations, Comm. Math. Phys. 214 (2000), no. 1, 137-189. 



18 



LAN-HSUAN HUANG, RICHARD SCHOEN, AND MU-TAO WANG 



[10] Corvino, J. and Schoen, R., On the asymptotics for the vacuum Einstein constraint equations, J. 

Diff. Geom. Volume 73, Number 2 (2006), 185-217. 
[11] Dain, S., Proof of the angular momentum-mass inequality for axisymmetric black holes, J. Diff. 

Gcom. 79 (2008), 3-67. 

[12] Huang, L.-H., On the center of mass of isolated physical systems with general asymptotics, Class. 

Quantum Grav. 26 (2009), no. 1, 015012, 25 pp. 
[13] Huang, L.-H., Solutions of special asymptotics to the Einstein constraint equations. Class. Quantum 

Grav. 27 (2010), no. 24, 245002, 10 pp. 
[14] Klainerman, S. and Nicolo, F., The evolution problem in general relativity. Progress in Mathematical 

Physics, 25. 

[15] Lindblad, H. and Rodnianski, I., Global existence for the Einstein vacuum equations in wave coor- 
dinates, Comm. Math. Phys. 256 (2005), no. 1, 43-110. 

[16] Regge, T. and Teitelboim, C, Role of Surface Integrals in the Hamiltonian Formulation of General 
Relativity, Ann. Phys. Volume 88 (1974), 286-318. 

[17] Zhang, X., Angular momentum and positive mass theorem. Comm. Math. Phys. 206 (1999), 137-155. 

Department of Mathematics, Columbia University, New York, NY 10027 
E-mail address: lhhuang@math.columbia.edu 

Department of Mathematics, Stanford University, Stanford, CA 94305 
E-mail address: schoenOmath.stanford.edu 

Department of Mathematics, Columbia University, New York, NY 10027 
E-mail address: mtwang@math.columbia.edu 



